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1.  INTRODUCTION 

Wyngaard  (1973)  introduced  the  concept  of  z-less 
stratification  for  cases  where  the  stratification  is  suffi¬ 
ciently  strong,  that  the  turbulence  no  longer  is  in  sig¬ 
nificant  communication  with  the  surface  (see  also  Holt- 
slag  and  Nieuwstadt,  1986).  Then  z  is  no  longer  a  pri¬ 
mary  scaling  variable,  nor  is  the  boundary-layer  depth. 
The  eddies  are  vertically  constrained  by  strong  stratifi¬ 
cation.  However,  the  z-less  concept  implies  more  than 
small  eddies,  since  vertically  continuous  turbulence  can 
still  organize  according  to  z  even  if  the  eddies  at  any 
level  are  small  compared  to  z-  For  example,  with  local 
similarity  where  the  relevant  Obukhov  length  must  be 
recast  in  terms  of  local  fluxes  at  level  z  instead  of  sur¬ 
face  fluxes  (Nieuwstadt,  1984),  the  overall  vertical  struc¬ 
ture  is  still  posed  in  terms  of  z/h  even  if  the  eddy  size  is 
small  compared  to  z.  In  this  sense,  local  similarity  still 
satisfies  the  criteria  for  traditional  boundary  layers.  On 
the  other  hand,  continuous  turbulence  between  the  sur¬ 
face  and  level  z  might  still  qualify  as  primarily  z-less  tur¬ 
bulence  if  the  principal  source  of  turbulence  is  detached 
from  the  surface  and  the  distance  above  the  ground  sur¬ 
face  is  only  a  secondary  influence. 

A  clear  example  of  z-less  turbulence  is  layers  of  tur¬ 
bulence  separated  by  nonturbulent  layer(s).  This  verti¬ 
cal  structure  does  not  satisfy  the  traditional  concept  of 
a  boundary  layer  in  that  surface-based  processes  are  of 
secondary  importance  (Mahrt,  1999).  Well-defined  labo¬ 
ratory  examples  of  such  “upside-down”  boundary  layers 
can  be  found  in  Ohya  (2(X)1).  In  this  study,  the  boundary 
layer  is  upside  down  if  turbulence  increases  with  height 
and  the  transport  of  turbulence  energy  is  downward  to¬ 
ward  the  surface.  While  the  upside-down  boundary  layer 
may  often  satisfy  the  conditions  for  z-less  turbulence, 
the  definitions  are  not  equivalent  in  that  the  upside-down 
class  requires  downward  transport  of  turbulence  energy. 

2.  DATA 

The  primary  data  set  is  6  levels  of  sonic  anemometer 
data  from  the  60-m  tower  in  CASES99  (Poulos  et  al., 
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2001;  Sun  et  al,  2002),  for  the  month  of  October,  1999 
over  relatively  flat  grassland.  In  addition,  we  include  data 
from  the  1.5-m  and  5-m  levels  of  a  mini  tower,  10  m  to 
the  side  of  the  main  tower.  The  sonic  anemometer  at  1.5 
m  was  moved  to  the  0.5-m  level  on  19  October  but  still 
included  in  the  study. 

The  data  was  quality  controlled  following  Vickers  and 
Mahrt  (1997).  Perturbations  are  defined  as  deviations 
from  a  stability-dependent  averaging  window  and  covari¬ 
ances  are  averaged  over  one  hour  to  reduce  flux  sampling 
errors.  Mostly  cloudy  cases  were  removed  by  discarding 
8  out  of  the  107  hourly  records  where  the  magnitude  of 
the  surface  net  radiation  loss  was  less  than  40  Wm^^. 

To  compute  mean  wind  profiles,  we  have  used  the 
R.M.  Young  propeller  anemometer  and  wind  vane  data 
for  15  m  and  above  and  the  sonic  anemometers  at  5  and 
1.5  m  (Figure  1).  The  profile  of  the  Richardson  num¬ 
ber  was  calculated  by  fitting  the  height  dependence  of 
wind  speed  and  potential  temperature  to  a  log-linear  pro¬ 
file  and  then  evaluating  the  gradients  analytically.  This 
Richardson  number  neglects  directional  shear,  which 
was  sometimes  important  near  the  surface. 

Vertical  transport  of  kinetic  energy  is  represented  by 
the  vertical  transport  of  vertical  velocity  variance.  The 
vertical  transport  of  the  horizontal  velocity  variances  be¬ 
haved  in  a  similar  manner  but  is  more  sensitive  to  the 
choice  of  averaging  time  and  to  random  flux  errors. 

3.  VERTICAL  STRUCTURES 

The  turbulence  energy  may:  1)  decrease  with  height 
across  the  tower  layer,  as  in  a  thin  traditional  shallow 
boundary  layer,  2)  vary  slowly  with  height,  as  in  a  bound¬ 
ary  layer  that  is  much  deeper  than  the  tower,  or  3)  in¬ 
crease  with  height  in  cases  where  the  turbulence  is  gen¬ 
erated  at  higher  levels  (upside-down  boundary  layer).  To 
separately  study  the  vertical  structure  of  these  states,  we 
composite  the  structure  of  the  turbulence  according  to 
classes  based  on  the  ratio  of  Cw  at  55  m  to  that  at  5  m. 
The  data  consists  of  107  one-hour  records  between  1900 
and  0600  local  standard  time  where  the  data  is  avail¬ 
able  for  all  eight  sonic  anemometers.  Eight  records  were 
eliminated  where  the  net  radiative  cooling  was  less  than 
40  presumably  due  to  clouds.  The  class  where 
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Figure  1:  The  composited  vertical  structure  of  the  a)  mean 
shear,  b)  local  gradient  Richardson  number,  c)  virtual  heat  flux 
and  d)  for  the  class  of  thin  traditional  boundary  layers 
(solid  lines)  and  the  upside-down  class  (broken  lines) 


where  the  ratio  exceeds  1.5,  which  captures  26  one-hour 
records.  The  thin  traditional  boundary  layer  where  the 
turbulence  decreases  with  height  will  be  defined  as  cases 
where  the  ratio  is  less  than  0.6,  which  captures  21  one- 
hour  records.  The  intermediate  class  corresponds  to  the 
case  where  the  turbulence  varies  slowly  with  height  and 
could  be  associated  with  a  boundary  layer  that  is  much 
deeper  than  the  tower.  For  the  intermediate  class,  the 
tower  is  not  tall  enough  to  capture  a  significant  frac¬ 
tion  of  the  boundary  layer  and  the  turbulence  quantities 
do  not  vary  substantially  across  the  tower  layer.  This 
class  also  includes  some  cases  of  poorly  defined  vertical 
structure.  This  intermediate  class  includes  52  one-hour 
records.  We  concentrate  on  the  vertical  structure  of  the 
thin  boundary  layer  and  the  upside-down  class  (Figure 
1). 

The  thin  traditional  boundary  layer  resolved  by  the 
tower  can  occur  any  time  during  the  night  but  most  of¬ 
ten  occurs  in  the  early  evening  with  very  stable  con¬ 
ditions  where  z/L  >  1  near  the  surface.  The  magni¬ 
tude  of  the  downward  heat  flux  in  the  composited  thin 
boundary  layer  decreases  linearly  with  height  to  approx¬ 
imately  30  m  and  then  is  small  above  this  level.  For  the 
thin  traditional  boundary  layer,  the  mean  shear  decreases 
monotonically  with  height  (Figure  la).  For  the  upside- 
down  class,  the  mean  shear  decreases  significantly  with 
height  only  up  to  20  m  and  increases  slightly  above  30 
m  becoming  twice  as  large  as  that  for  the  thin  traditional 
boundary  layer  at  the  40-m  level.  For  individual  cases, 
the  local  shear  in  the  upper  part  of  the  tower  layer  may  be 


locally  large  but  the  level  of  the  significant  shear  varies 
between  cases,  and  does  not  occur  at  any  level  for  some 
cases.  Therefore,  the  composited  shear  profile  does  not 
show  any  well-defined  shear  maximum. 

The  greater  shear  for  the  upside-down  class  at  higher 
levels  exerts  a  strong  influence  on  the  Richardson  num¬ 
ber  at  these  levels,  partly  because  the  Richardson  number 
depends  quadratically  on  the  shear  and  partly  because  the 
stratification  decreases  rapidly  with  height  (not  shown). 
The  Richardson  number  increases  with  height  within  and 
above  the  thin  traditional  boundary  layer  (Figure  lb), 
although  the  Richardson  number  profile  is  sensitive  to 
the  method  of  calculation.  For  the  upside-down  class, 
the  Richardson  number  does  not  increase  with  height 
above  30  m,  in  contrast  to  the  traditional  thin  bound¬ 
ary  layer.  While  the  buoyancy  flux  decreases  to  immea¬ 
surably  small  values  above  30  m  for  the  thin  traditional 
boundary  layer,  the  magnitude  of  the  buoyancy  flux  in¬ 
creases  above  30  m  for  the  upside-down  class.  The  ver¬ 
tical  structure  of  the  stress  (not  shown)  is  similar  to  that 
of  the  buoyancy  flux  but  slightly  less  organized. 

For  the  thin  traditional  boundary  layer,  the  vertical 
transport  of  vertical  velocity  variance  is  near  zero  or 
weak  upward  (Figure  Id).  Since  the  length  scale  of  the 
turbulence  is  small  for  strongly  stratified  boundary  layers 
with  weak  turbulence,  the  turbulence  is  more  efficiently 
dissipated  locally.  The  downward  transport  of  the  verti¬ 
cal  velocity  variance  for  the  upside-down  class  (Figure 
Id)  implies  that  this  class  is  not  simply  a  boundary  layer 
which  is  much  deeper  than  the  tower  layer,  but  contains 
turbulence  associated  with  shear-generation  that  is  not 
directly  coupled  to  the  surface.  This  is  not  a  traditional 
boundary  layer. 

4.  MIXING  LENGTH 

Because  boundary-layer  depths  are  often  20  m  or  less, 
the  surface-layer  prediction  of  the  mixing  length,  K  z/(!>m 
is  often  expected  to  be  valid  at  night  only  in  the  low¬ 
est  few  metres.  The  failure  of  similarity  theory  at  higher 
levels  is  supported  by  the  composited  diurnal  variation 
of  the  aerodynamic  roughness  length  (Figure  2).  The 
roughness  length  was  computed  using  Monin-Obukhov 
similarity  theory  and  the  observed  heat  and  momentum 
fluxes.  At  1.5  m,  the  roughness  length  is  approximately 
independent  of  time  of  day.  However,  for  the  5-  and  10- 
m  levels,  the  aerodynamic  roughness  length  dramatically 
decreases  at  night.  Since  the  physical  roughness  of  the 
surface  does  not  change,  this  variation  is  due  to  inappli¬ 
cability  of  surface  layer  similarity  theory.  As  a  conse¬ 
quence,  the  study  of  the  surface  layer  similarity  theory 
requires  eddy  correlation  measurements  within  the  low¬ 
est  few  metres.  Note  that  the  thermal  roughness  length 
based  on  the  surface  radiation  temperature  varies  sub- 


Figure  2:  The  dependence  of  the  composited  roughness  length 
on  time  of  day  for  the  1.5>  5-  and  10-m  levels. 


stantially  even  at  1.5  m  over  this  simple  surface. 

However,  formulating  the  mixing  length  in  terms  of 
based  on  local  fluxes  at  level  z  provides  a  credi¬ 
ble  approximation  to  the  mixing  length  even  above  the 
surface  layer  up  to  about  40  m.  Above  40  m,  this  ap¬ 
proximation  overestimates  the  mixing  length.  Part  of  the 
success  above  the  surface  is  that  this  approximation  be¬ 
comes  approximately  z-less  for  strong  stability  (Mahrt, 
1999).  However,  part  of  the  success  is  due  to  self  corre¬ 
lation  in  that  both  the  mixing  length  and  Kz/(t>m  are  pro¬ 
portional  to  M». 

Dismissing  similarity  theory  for  all  but  the  lowest 
tower  level,  we  now  explore  the  relationship  between  the 
mixing  length  and  the  gradient  Richardson  number  (Fig¬ 
ure  3).  Here  we  focus  on  the  5-m  level,  which  is  almost 
always  within  the  boundary  layer,  and  the  40-m  level, 
which  is  generally  above  the  boundary  layer  except  dur¬ 
ing  strong  mixing  periods.  Several  features  transcend 
the  behavior  of  the  mixing  length  at  all  of  the  levels. 
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Figure  3:  The  dependence  of  the  mixing  length  on  the  layer 
Richardson  number  centered  about  the  eddy  correlation  level 
and  the  dependence  of  the  mixing  length  on  the  bulk  length 
scale  (Eq.  2)  for  the  5-  and  40-m  levels. 


whether  it  is  within  the  boundary  layer  or  above  it.  The 
mixing  length  decreases  sharply  from  near-neutral  condi¬ 
tions  (Ri  0)  to  modestly  stable  conditions  (Ri  w  0.20). 
The  relationship  is  rather  successful  considering  large 
random  flux  errors  and  difficulties  of  estimating  mean 
vertical  gradients. 

For  stronger  stability  (Ri  >  0.20),  the  value  of  the 
mixing  length  tends  to  remain  constant  with  large  scat¬ 
ter.  Surprisingly,  eliminating  cases  with  large  intermit- 
tency,  nonstationarity  and  downward  transport  of  turbu¬ 
lence  energy  does  not  visibly  reduce  the  scatter.  The 
intermittency  is  smallest  near  the  surface  and  increases 
with  height. 

Near-neutral  conditions  include  numerous  cases  of 
significant  nonstationarity.  These  cases  are  associated 
with  transition  periods  where  the  flow  changes  between 
unstable  and  stable  or  between  stable  and  unstable. 
Nonetheless,  the  value  of  the  observed  mixing  length 
at  neutral  conditions  is  consistent  with  the  asymptotic 
trends  of  the  stability  dependence  for  the  stable  and  un¬ 
stable  regimes.  The  near-neutral  value  of  the  mixing 
length  increases  with  height  up  to  about  20  m  and  then 
becomes  approximately  constant  with  height  with  a  value 
of  6  m.  The  boundary-layer  depth  does  not  influence 
the  height  dependence  of  the  neutral  mixing  length  for 
the  composited  profiles  in  an  obvious  way  possibly  be¬ 
cause  the  boundary-layer  depth  varies  substantially  be¬ 
tween  the  different  records  and  is  not  well  defined  for 
many  cases. 

To  model  the  z-less  mixing  length,  we  pursue  a  simple 


dependence  on  the  Richardson  number 

=  lo{z)exp(~-aRi)  (1) 

The  principal  closure  problem  is  specification  of  the  mix¬ 
ing  length  at  neutral  stability,  Iq.  For  the  present  data, 
to  =  6w  is  an  adequate  approximation  but  the  general¬ 
ity  of  this  value  is  not  known  and  the  specification  of  a 
dimensional  quantity  is  unappealing. 


4 A  Alternative  formulation 


The  principal  disadvantage  of  Eq.  1  is  specification  of 
the  asymptotic  value  at  neutral  stability.  This  difficulty 
can  be  avoided  by  formulating  the  mixing  length  as 


S{z)^C(z) 


{^/e)A© 


(2) 


where  C{z)  is  a  nondimensional  coefficient.  This  rela¬ 
tionship  explains  a  substantial  fraction  of  the  variance  of 
the  mixing  length  and  contains  no  self  correlation.  The 
coefficient  C{z)  increases  with  height  and  asymptotes  to 
a  roughly  height-independent  value  of  2.5  x  10*"^.  The 
disadvantage  of  this  approach  is  that  it  asymptotes  to  in¬ 
finity  at  neutral  conditions  where  A©  vanishes,  at  least 
for  a  fixed  Az.  However,  with  increasing  mixing  length, 
the  relevant  temperature  difference  needs  to  be  com¬ 
puted  over  larger  scales  corresponding  to  the  eddy  size, 
in  such  a  manner  that  A©  recognizes  the  stratification  at 
the  edges  of  the  mixing  layer,  or,  becomes  constrained 
by  the  ground  surface. 

During  the  oral  presentation,  several  different  models 
which  combine  surface  layer  similarity  theory  and  the 
above  z-less  relationship  will  be  evaluated. 
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